Thermal memory: a storage of phononic information 



OO 
O 

o 

(N 

ojo: 
< 

in 



I 

C 

o 
o 



o 

OO 

o 



X 



Lei Wang 1 ' 2 and Baowen Li 2,3 * 
1 Department of Physics, Renmin University of China, Beijing 100872, P. R. China 
2 Department of Physics and Centre for Computational Science and Engineering, 
National University of Singapore, Singapore 117542 
3 NUS Graduate School for Integrative Sciences and Engineering, Singapore 117597, Republic of Singapore 

(Dated: 25 August 2008) 

Memory is an indispensable element for computer besides logic gates. In this Letter we report 
a model of thermal memory. We demonstrate via numerical simulation that thermal (phononic) 
information stored in the memory can be retained for a long time without being lost and more 
importantly can be read out without being destroyed. The possibility of experimental realization is 
also discussed. 
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Unsatiable demands for thermal management/control 
in our daily life ranges from thermal isolating to efficient 
heat dissipation have driven us to better understand heat 
conduction from molecular level. Recent years has wit- 
nessed a fast development including theoretical proposals 
in functional thermal devices such as thermal diode that 
rectifies heat current [l| , thermal transistor that switches 
and modulates heat current 0, heat pump that car- 
ries heat against temperature bias and experimental 
works such as nanotube phonon waveguide [J], thermal 
conductance tuning [|| and nano-scale solid state ther- 
mal rectifier |6[ . More importantly, logic operations with 
phonons/heat have been demonstrated theoretically 0, 
which, in principle, has opened the door for a brand new 
subject - phononics - a science and engineering of process- 
ing information with heatjij]. The question arises natu- 
rally and immediately: whether a thermal (phononic) 
memory that can store information is possible? 

In this Letter, we would like to give a definite an- 
swer by studying the transient process, which in fact 
exhibits much richer phenomena than asymptotical sta- 
tionary state does, of a non-equilibrium system. This 
topic has been however rarely studied so far. 

Like an electronic memory that records data by main- 
taining voltage in a capacitance, a thermal memory stores 
data by keeping temperature somewhere. An ideal mem- 
ory keeps the data forever without fading. This is never 
possible for a thermal system because sooner or later 
randomness of the heat (fluctuation) will eliminate any 
record of the history. However this problem is not se- 
rious from application point of view since we only need 
to maintain the data until we refresh it or read it. This 
is exactly the case for the widely used Dynamic Ran- 
dom Access Memory (DRAM) that needs refresh opera- 
tion regularly. Any system that keeps temperature (thus 
data) somewhere for a very long time might be a candi- 
date for thermal memory, such as breather and localized 



harmonic mode etc. However perturbation is unavoid- 
able in a thermal system, especially as the data are read 
(i.e., the local temperature is measured). The breather 
and harmonic mode, even if do not collapse immediately 
under the perturbation are not able to recover their 
original state. Because the energy that is changed by 
the perturbation, e.g., the necessary energy exchange in 
order to build equilibrium between the system and the 
thermometer during the data reading (temperature mea- 
suring) process, is not recoverable in the autonomous sys- 
tems without an external energy source/sink. We thus 
have to turn to a thermal-circuit with power supply, i.e., 
driven by external heat bath. 

Like an electric circuit, the temperature and heat cur- 
rent distributions of a thermal-circuit are determined by 
Kirchhoff's laws. In any linear thermal-circuit, i.e. all 
thermal resistances are fixed, independent of tempera- 
ture and/or temperature drop, the steady state that sat- 
isfies Kirchhoff's laws must be unique. This unique state 
must be stable, namely under any perturbation the sys- 
tem will eventually return to this state. However in or- 
der to work as a thermal memory, a thermal-circuit must 
have more than one stable steady states, such a bi-stable 
thermal-circuit is only possible in the presence of nonlin- 
ear thermal devices. 

To make such a thermal-circuit, we consider a one di- 
mensional nonlinear lattice that consists of two Frenkel- 
Kontorova (FK) segments [l(| sandwich a central particle 
in between. Its configuration is shown in FigfTJa) and its 
Hamiltonian follows: 



H = H L + H h 
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and the Hamiltonian of each segment L/R can be writ- 

i=l — 2^- + —£ i! -{XL/R,i-XL/R,i-l) ~ 

(27rp uua z,7rx L/R,ij with x^/ft i denotes the displacement 
from equilibrium position of the i th particle in segment 
L/R. We have set the mass of each particle to unity. 
The parameters k and V are the harmonic spring con- 
stant and the on-site potential of the FK segments. 
We couple the last particles of segments L and R to 
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particle O via harmonic springs. Thus Hi n t = -§- 



T,l/r M 2 >R ( x l/r.n l/r ~ x ) 2 ~ COS27TX . In our 
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numerical simulation, segment L and R each contains 50 
particles. Fixed boundary conditions are applied. Other 
parameters are: K L =K intL =0A5, K R =0.2, K int R=0.05, 
Vl=10, Vr=0, Vb=10. The deep valley of on-site po- 
tential of particle O, Vo, and the weak coupling between 
segment R and particle O, Ki n tR, induce a negative dif- 
ferential thermal resistance (NDTR)Q] between O and 
segment R, which is the key for the model of thermal 
memory. We use a wavy curve to stress this in FigJIJa). 
Tl and T#, temperatures of the power supplies (heat 
baths) that contacted to the ends of the two segments, 
are fixed to 0.04 and 0.3, respectively. Without these en- 
ergy resource/sink, the system reduces to an autonomous 
one thus, as has been explained above, is not able to work 
as a thermal memory. All heat baths are simulated by 
Langevin heat baths 111 ] . 

Because of NDTR, with fixed Tl/Tr and a control 
heat bath with adjustable temperature To is coupled to 
the particle O, there exist more than one values of To 
with which the heat current from control heat bath to 
particle O, Jo, is zeroQ, as shown in FigQJb). Suppose 
the temperature of O is initially set to either "T on " or 
"T ff" by the control heat bath. After the control heat 
bath is removed the state will remain unchanged for a 
long time, in spite of the thermal fluctuation. Because as 
the consequence of the small perturbation which slightly 
changes the temperature To, the changes of heat currents 
5Jl and SJr always pull To back, thus these two states 
are "stable" . It is easy to check the other steady state in 
between is however unstable. 

In the equilibrium state of an Hamiltonian system with 
crgodicity, temperature can be defined as twice of the 
average kinetic energy per degree of freedom. In or- 
der to study transient process, the "finite time temper- 
ature" is naturally defined as, for a ID system, twice 

of the average kinetic energy in a time window St, i.e.,: 
t+ — 

T(t) = / t _a? v 2 (t')dt' . Hereafter, any temperature 

which is a function of time t means the finite time temper- 
ature. Due to the thermal fluctuation, even in a steady 
state, finite time temperature has different value at dif- 
ferent time window. The probability density function 
(PDF) of the finite time temperature is thus defined ac- 
cordingly. If the time window St is chosen to be too 
long, the PDF approaches a S function locates at "infi- 
nite time" temperature T, while if St is too short, big 
thermal fluctuation dominates thus the distribution is 
more or less a Maxwell energy distribution. In the me- 
dial cases PDF presents much more information other 
than the value of temperature T itself. In this Letter, St 
is fixed to 10 4 dimensionless time units. Since the fre- 
quency of the system we studied is about 0.2 this 
period of time covers about 2,000 oscillations, thus is a 
microscopical long time. 



Since "on" and "off" are two stable states, the finite 
time temperature of the particle O in the absence of con- 
trol heat bath is expected to stay around T on and T Q ff 
with much higher probability than elsewhere. In FigQJc) 
the locations of the two peaks clearly confirms this expec- 
tation. By adjusting parameters of the system, the values 
of T on and T Q f f and the peaks shift simultaneously but 
always coincide with each other. 
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FIG. 1: (Color-online) (a) Configuration of a thermal mem- 
ory. Black wavy curve indicates the weak coupling between 
particle O and segment R, which results in the negative dif- 
ferential thermal resistance (NDTR) here, (b) Heat currents, 
Jl, Jr and Jo versus temperature To- The NDTR between 
particle O and segment R makes Jr increases in a wide region 
when To is increased, thus makes the curve cross with Jl at 
more than one values of temperature To ■ (c) PDF of the fi- 
nite time temperature To(t) of the central particle O in the 
absence of control heat bath. Notice the clear correspondence 
of the two stable steady states in (b) and the two peaks in (c), 
which are indicated by the two vertical dot lines. Since To (i) 
can only change continuously, the very low PDF in between 
clearly points out the low transition probability, thus implies 
the stability of the two states "on" and "off" . 

We have so far confirmed the stability of the two states 
"on" and "off" of the model. In the following, we shall 
demonstrate a complete, four-stage, writing-reading pro- 
cess of the thermal memory, as is shown in Fig(5Ja) . Dur- 
ing the whole process the left and right ends of the mem- 
ory are always connected to two power supply heat baths 
with temperatures Tl—0.04 and Tr=0.3, respectively. 
Stage 1: initializing the memory. Each particle is 
initially set a velocity chosen from a Gaussian distribu- 
tion corresponding to a temperature in (T l ,Tr). After 
some time the system approaches to an asymptotic state, 
thus To(t) saturates, see FigJ^b). Stage 2: writing 
the datum into the memory. In our simulation, the 
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"writer" is simulated by a FK lattice with 7V=10 particles 
and identical parameters with segment L. It is connected 
via a linear spring with fc=1.0, at its one end, to the par- 
ticle O of the memory. The temperature of the writer is 
initially set to T f / and the other end is connected to a 
heat bath with temperature T ff also. In a short time 
the writer cools down To(t) to T a ff, namely the datum 
"off" has been written into the memory. Stage 3: main- 
taining the datum in the memory. In this stage, the 
writer is removed from the memory. And we can see 
from Fig(2Jb) that To{t) remains nearly unchanged dur- 
ing this stage, which means that the datum stored in the 
memory can last for a long time. Stage 4: reading out 
the datum from memory. A "reader" (thermometer) 
is used to read the datum out from the memory. Here 
the "reader" is simulated by the same FK lattice as the 
writer. It is connected, to the particle O of the mem- 
ory. The temperature of the reader is initially set to a 
middle temperature between T on and T ff ( m this Let- 
ter it is 0.11). Different from the writer, the reader is 
not connected to any heat bath. This stage is the most 
important and most interesting one. We see that at the 
beginning the particle O is temporarily warmed up be- 
cause the reader is hotter. However, since the "off" state 
is stable, To(t) and also the temperature of the reader 
are cooled down to T a f f shortly, implying that the datum 
has been successfully read out. To gain a clearer picture, 
we show the heat currents through the two segments Jl 
and Jr in the Figj2fc). At the beginning of the reading 
stage since the temperature of particle O is warmed up 
by the reader, as a response of the increase of To(t), Jl 
increases a lot whereas Jr increases only very little, thus 
the net current Jr-Jl changes from zero to negative, i.e., 
the power supply (left end heat bath) absorbs heat from 
the particle O. This is the reason that To can be cooled 
down and recovered to T Q f f automatically. 

The results shown in Figj2] are obtained from an en- 
semble average over 20, 000 independent samples. For a 
model of thermal memory, the ratio of samples that keep 
the correct data is an important feature. To describe 
this quantitatively, we calculate the ratio of "on" state: 
r on (t) = N on (t)/N, while N on (t) is the total number of 
samples stay at "on" state at time t and N is the total 
number of samples. We define, at time t, a sample is 
at "on" state if the finite time temperature of its parti- 
cle O is greater than a critical value T c =0.11. we show 
r n{t) during the whole process in FigOfd). It is clearly 
seen that after the writing stage, r on (t) always keeps at 
a very low level, even under the not-very-small pertur- 
bation from the hot reader. The discrepancy between 
r on (t) and is in the order of 10~ 4 , which means two or 
three "errors" among the 20,000 samples, is even indis- 
tinguishable in the figure by eye, namely the memories 
successfully maintain the data "off" that are initially set 
by the writer, the data are hardly destroyed even after 
being read! 
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FIG. 2: (Color-online) A writing-reading process for an 
"off" state. Vertical dot lines separate different stages to 
guide eyes. (a) Illustration of the four stages, (b) Tempera- 
tures of the central particle To(t), writer T mr i t er(t) and reader 
T rea der{t)- (c) Heat currents through the two segments, Jl 
and Jr. It can be seen that at the beginning of the reading 
stage, the net current Jr-Jl changes to negative, namely, en- 
ergy is absorbed from O thus To(t) is cooled down back to 
the low temperature T //. (d) After the writing stage r on (t) 
keeps very close to zero, namely nearly all samples stay at the 
"off" state, even under the not-very-small perturbation from 
the hot reader. 



It has been pointed out that, to be a memory, the sys- 
tem must have more than one stable steady states around 
which the above writing-reading process can be com- 
pleted, thus in the following we demonstrate the writing- 
reading process around the other stable steady state, the 
"on" state in Fig|3] The difference from FigJ5] is that 
in the second stage the writer is set to temperature T on 
and connected to a heat bath with temperature T on . We 
see in this stage To{t) changes to T on shortly. In the 
consequent maintaining stage To(t) keeps unchanged at 
T on . In the reading stage: the reader is initially set to 
the same temperature as before. We see at the begin- 
ning because this time the reader is colder, the particle 
O is temporarily cooled down, as a response Jl decreases 
while Jr increases, thus the net current Jr-Jl becomes 
positive, which warms up the particle O (also the reader) 
and recovers its temperature to T on . The datum "on" is 
again read out successfully. As for r on (t), see Fig[3£d). 
This time r on (t) keeps at a very high level. The discrep- 
ancy between r on (t) and 1 is again at the order of 10~ 4 , 
namely the datum "on" can be successfully maintained 
in and read out from memory. 
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FIG. 3: (Color-online) A writing-reading process for an 
"on" state, (a) Illustration of the four stages, (b) Tempera- 
tures of the central particle To(t), writer T mr it er (t) and reader 
T rea der{t)- (c) Heat currents through the two segments, Jl 
and Jr. At the beginning of the reading stage, the net current 
Jr-Jl is positive, namely provides energy to O thus warms 
it up back to the high level T on . (d) After the writing stage 
Tonit) remains very close to one, namely, nearly all samples 
stay at the "on" state, even under the not-very-small pertur- 
bation from the cold reader. 



variance a which describes fluctuation of the heat cur- 
rent. The mean first passage time (MFPT) through a 
certain boundary X = AX (which denotes T c — T a ff 
or T on — T c ) conditional upon X(0) = represents the 
average lifetime of the data in the memory. In the lim- 
iting case that MFPToc^erfi(^i), where 
etS.(x)=J^ e x dx' [lj]. MFPT diverges rapidly as a 2 



decreases. Parallel combining n identical memories to- 
gether decreases a by n times while keeps A unchanged, 
thus increases the average lifetime of data greatly. Such 
a fast divergence tendency of MFPT has also been widely 
found in a class of diffusion processes with steady state 
distribution [1 51 ] , thus the above conclusion is not limited 
to the specified process, the OU process. 

In summary, we have demonstrated the feasibility of 
thermal memory from a bistable thermal-circuit. The in- 
formation stored in such a thermal memory can last very 
long time and, more importantly it is self-recoverable un- 
der the not-very-small perturbation from the thermome- 
ter when the data are read. With the rapid developing 
nano-technology @, H, @] , the thermal memory should be 
realized in nanoscale systems experimentally in a foresee- 
able future. 
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000-203-112 from Ministry of Education of Republic of 
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Finally we would like to discuss the lifetime of data in 
the thermal memory. In an electric DRAM, since real 
capacitors leak charge, the data will eventually fade un- 
less the capacitor charge is refreshed periodically (about 
~ 100/s). Similarly, we find a few "errors" among the 
many samples that we studied at the end of the writing- 
reading process. Suppose the data fading process is a 
Poisson process then the average lifetime of the data is 
roughly 5 * 10 9 dimensionless time units which contains 
about 10 9 periods of oscillation. This corresponds to 
100 /is if the memory is made of carbon nanostructures 
such as nanotube[l2j, in which the key of the model, 
NDTR has been found numerically [liij. This lifetime is 
unsatisfying when compared with an electronic DRAM. 
However this value can be easily and greatly enlarged 
by parallel combining more identical memories together. 
The dynamics of finite time temperature of particle O, 
To (t) around either of the two stable steady states can be 
roughly described by an autoregressive diffusion process, 
say, the simplest one, Ornstein-Uhlenbeck (OU) process: 
X = -XX + er£(i). Where X corresponds to T (t) rel- 
ative to T on or T ff. A>0 represents the effect of the 
response net heat current that pulls X back to zero and 
cr£(i) is a Gaussian white noise with zero mean and fixed 
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